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Spontaneous Pair Creation revisited 
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Recently the so called Spontaneous Pair Creation of electron positron pairs in a strong external 
field has been rigorously established. We give here the heuristic core of the proof, since the results 
differ from those given in earlier works. 
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I. INTRODUCTION 



The creation of an electron positron pair in an almost 
stationary very strong external electromagnetic field 
has been called — unfortunately misleading — spontaneous 
pair creation (SPC) (0 - The phenomenon 

emerges straight forwardly from the Dirac sea interpreta- 
tion of negative energy states: An adiabatically increas- 
ing field lifts a particle from the sea to the positive en- 
ergy subspace where it scatters and when the potential 
is gently switched off one has one free electron and one 
unoccupied state — a hole — in the sea. The experimen- 
tal verification has been sought in heavy ion collisions, 
but without success so far A rigorous proof 

of the existence of SPC has been lacking until recently 
[T^ |. Surprisingly our rigorous proof — the heuristics of 
which we give here — yields time estimates different from 
what has been reported earlier and which may be rele- 
vant for experiments. The rigorous result concerns the 
so called external field problem, i.e. interactions between 
the charges are neglected. Vacuum polarization will in 
general perturb the external field (see 01 f° r a rigorous 
attempt) and one may think of the external field as an ef- 
fective field. The description of SPC in second quantized 
external field Dirac theory is known to be equivalent to 
the existence of certain type of solutions of the one par- 
ticle Dirac equation. We shall only discuss the latter. 



II. FORMULATION OF THE PROBLEM 



Consider the one particle Dirac equation with external 
electric potential A, a real valued multiple of the 4x4 
unit matrix. Amc 2 gives the potential in the units eV. 



On microscopic timescales r = ls f-t the equation reads 



.dip T 



■ ^2 a l°l^r + A eT (x)i' T + /3Vr 



1=1 



(1) 



where s is a parameter representing the slow time varia- 
tion of the external potential. 

We introduce in (JIJ the macroscopic time scale s = er: 



os e 



(2) 
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We will describe pair creation using the Dirac sea in- 
terpretation of the Dirac equation. The spectrum of 
the Dirac equation without external field is (-co, — 1] n 
[1, oo). Wave functions which lie in the corresponding 
positive energy subspace of the free Dirac operator are 
interpreted as wave functions of particles or electrons. 
The Dirac sea is a many particle wave function built out 
of wave functions of negative energy (states) of the free 
Dirac operator. In the so called vacuum all states in the 
sea are occupied by particles. "Holes" in the Dirac sea are 
unoccupied states which are interpreted as anti-particles 
or positrons. 

Suppose now that the potential is adiabatically 
switched on and later switched off. When the potential 
is non zero there may be bound states in the gap, [— 1, 1]. 
The eigenvalues of these bound states change with the 
strength of the potential. The adiabatic theorem (see 
e.g. ensures that there is no tunnelling across spec- 
tral gaps. So as long as all bound states are isolated 
^from the upper continuous spectrum, transitions from 
the negative to the positive energy subspace are adiabat- 
ically not possible and the probability of creating a pair 
is zero (undercritical case). When the external field be- 
comes overcritical (at time s c — 0), i.e. when a bound 
state vanishes in the upper energy subspace, there exists 
a solution of the Dirac equation which follows the path of 
this bound state from the negative to the positive energy 
subspace. If the wave function (or a part of it) does not 
"follow" the bound state back into the negative energy 
subspace when the potential is switched off again, but 
stays in the positive energy subspace, there will be a hole 
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in the sea and a freely moving particle: SPC is achieved. 
The scenario is symmetric under change of sign of the po- 
tential: It then transports an unoccupied state (a hole) 
from the positive energy subspace to the sea and catches 
a particle from the sea when it is switched off. The hole 
(positron) then scatters. 

To show to what extend the scenario in fact holds one 
must control the propagation of the wave function emerg- 
ing from the bound state during over-criticality. The 
wave function will generically not be a nice scattering 
state, i.e. its momentum distribution will be large for 
small momenta, that is, the wave function lingers around 
the range of the potential rather than moving away. We 
shall describe now the heuristic core of the situation when 
SPC happens. 



with imaginary k — in. In the following we will assume 
that A is such that for each time a there exist either no 
or one (degenerate) bound state and that the bound 
state vanishes in the positive continuum at a = 0. Since 
there generically exists a bound state <E> with energy 1 
(i.e. k = 0) at time a = and d a E a \ a= o ^ (see [l3| 
and our argument below), we will only discuss that case. 

Using for a — k — observing the explicit form of 
Gq (see e.g. 0) one can easily estimate $ for large x. 
Assuming $ bounded, a term falling off like x~ l appears. 
Using that this term has to be equal to zero (so that 
is square integrable) one gets the crucial identity 

J(l + f3)A (x)<S>( X )d 3 x = . (8) 



III. THE SPECTRUM OF THE DIRAC 
OPERATOR 

We expand the wave function in generalized eigenfunc- 
tions which by themselves are time dependent. We shall 
need the eigenfunctions for times a close to the critical 
time. Consider the eigenvalue equation 

Ecp = D a (p (3) 

for fixed <r€l. The continuous subspace determined by 
k £ K 3 \{0} or k = and cr ^ is spanned by generalized 
eigenfunctions (not square integrable) ip J (k, cr, x) , j — 
1,2,3,4, with energy E = ±E k = ±y/k 2 + 1, j = 1,2 
being solutions with positive energy. For ease of notation 
we will drop the spin index j in what follows. 

The generalized eigenfunctions also solve the Lipp- 
mann Schwingcr equation 

<^(cr,k,x) = (y3 (k,x) 

+ J G+(x - x')A a (x'Ma, k, x')dV , (4) 

with ipo(k, x) = £(k)e* I! K £k ' x , the generalized eigenfunc- 
tions of the free Dirac operator D°, i.e. G k is the kernel 
of (E k - D )- 1 = lim^o(£fc -D° + iS)- 1 0. 
Introducing the operator 

T k a j = J G+(x - x').4 CT (x')/(x')<iV, (5) 

(0} becomes 

(l-T ff fe M<T,k,-)=^o(k,-)- (6) 

There may also exist resonance-eigenfunctions which are 
not square integrable but which decay as x — > oo. 

Finally there may exist eigenfunctions, i.e. bound 
states $ CT (re) of © with energies E = y/l — k 2 e [—1,1], 
They satisfy instead of iJBJ 

(1-T CT K (K) )$ CT(K) =0, (7) 



A. Propagation estimate 

We estimate the propagation of a wave function gener- 
ated by the static Dirac Operator D a = D° + A a (x), 
where a > should be thought of as near the criti- 
cal value (the relevant regime turns out to be of order 
< 7 = 0(e 1 /3)). 

Since the generalized eigenfunctions for (cr, k) rj (0,0) 
are close to the bound state $ it is reasonable to write 
in leading order: 

^(a,k,x)«77 a (k)$(x) . (9) 

Since they solve (0J, the first summand of (0J must be- 
come negligible with respect to 77 a (k)$, which is part 
of the second summand. Hence ^(k) must diverge for 
(cr, k) — > (0,0). For the outgoing asymptote of the state 
$ (generalized Fourier transform) evolved with D a near 
criticality we have with @ that 

$ out (er,k) := J (2ir)-i<P(x)lp(a,k,x)d3x 

« (27r)-i^(k). (10) 

Now, for (a, k) close to but different from (0, 0), ^(k) ~ 

$ out (cr, k) will be peaked around a value k(a) with width 
A(cr) (determined below) defined by 

Va (k(a)±A(a))^r ]rT (k(a))/V2 . (11) 

We may use the width for the rough estimate 

\d k $ out (a,k)\ <A(cr)- 1 $ 0ttt (o-,A tr ), (12) 

where the right hand side should be multiplied 
by some appropriate constant which we — since it is 
not substantial — take to be one. Using JJjJ, QlOp. 
d 3 k — k 2 dfldk and partial integration (observing 

^dfce^+^f = e-^+^i) we get 

C/ CT ( S ,0)$ = e - isD "$ m — ^ [ e-^+^i^outipdSk 
= ^ J e-« 1+ ^id k (\$ oat \H(x)kdn) dk 
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By HU, and JTTJ, 

assuming that A(cr) fc(c) 



Expanding and #3 around a = we thus obtain 



A(cr)fc fc 2 

-d 3 fc 



A(cr)fc(cr) 



Hence 



|f7 CT (s,0)$(x)| < 



2g[$(x)l 
sA(ct)A;(ct) 



|$ OMt | 2 d 3 fc. 



Since 4> is normalized we get for the decay time Sd, de- 
fined by \{U(s d ,0)$,®)\ « 1/2 



Sd « 4e(fc(a)A(a))- 



(13) 



B. Control of the Generalized Eigenfunctions 

Let us now estimate r] a (k) for u « 0. In view of © 
and © we have that 



(1-T^> ff (k)$«^ (k,-) 



(14) 



We can estimate r\ a (k) by considering the scalar product 
of witn ^0$: 

»j a (k)((l- 1^)*,^,*) « ^ (k,-),A $) . 

Since yo(0, x) = £(0) for which 

^(0) = £°£(0) = /3£(0) (15) 

we have that y>o(0,x) = ^p£(0) and in view of (JSJ 

fco, A,*) |fc= = ffl- /" (1 + /?)A (x)$(x)d 3 z = . 

(16) 

Hence (tpo(\t, •) , Ao$)) = O(k), and it turns out that 
0(k) = Ck + 0{k 2 ) with an appropriate C ^ 0. Thus 

r^^Ckdl-T^A^))- 1 

Expanding T% in orders of k around k = until fourth 
order yields 

r] a (k) m Cfc (S£ + fcSJ + fc 2 S 2 °" + fc 3 ^ + **S£) _1 .(17) 

By virtue of JSJ one easily shows 

(T>, = (18) 

and we get with J7J, using A CT > Aq, that 

S£ = ((l-T 17 0)*,A)*) = <(TOo-T ff O)$,i4o$) 
= ((A -A a )<!>,T 03>) = {(A -A ( ,)<i>,$) 

= - || ^jA a - A $ || 2 = -C7 ff • 

with Co > 0. Computing that dkG^ \k=o= 1 + P, and 
observing (JSJ and ijHJl, duT k <& \k=o= 0, hence with {THJ) 

5f - %(T*$,i4o*) |*=o= 9 fc (T fc $, | fc=0 = . 



»7«r(k) 



(19) 



-Ck (C <7 + (C 2 + ©((t))^ 2 + (C* 3 + G(ct))A; 3 



We shall now determine the constants C 2 and C3. For 
that we do a similar expansion for J7J, i.e. for a < 
and for fc = j/t. Dividing Ijl4(l by rj a {k.) and replacing $ 
by $<t( k )(~ 4) yields the left hand side of and thus 
'7<t(k)(* k ) = °o (otherwise the right hand side of J7J would 
not be zero). Therefore for all k 



-C \a(K)\~(C2+O(a(K)))K 2 ~i{C 3 +O(a(K)))K 3 



= 



Since C > 0, assuming C 2 , C 3 ^ (see jT^I for a proof) 
we conclude that C 2 < and C3 must be imaginary. 
Note that near k = 0, a ^ k 2 = E — m, while if C\ 7^ 0, 
E — m ~ — a 2 , i.e. dE/da = at cr = 0, suggesting that 
a transversal crossing of the eigenvalue into the upper 
continuum goes together with a bound state at £7 = 1. 
Hence for l|19(l we get 



-Cfc 



Cod - (|C 2 | + 0{o))k 2 - i(\C 3 \ + 0(a))fc 3 



For Cocr « C 2 fc 2 the denominator behaves like 6*3/03, oth- 
erwise it behaves like Cqg — C2k 2 . Hence by (|1(J|) 

\$ out (<r,k)\ 2 « cfc 2 ((Coct - |c 2 |fc 2 ) 2 + |c 3 | 2 fc 6 ) 1 (20) 

This result 0| differs from the results given in the liter- 
ature (see e.g. formula (7) in Q). The right hand side 
of (|2U[I obviously diverges for (<r, fc) — > (0,0). For fixed 
0^(T~0 the divergent behavior is strongest close to the 
resonance at (Cqct — |C 2 |fc(er) 2 = 0) 



k{a) = y/aColCil' 1 = 0(VF) . (21) 

In view of (| 1 1 11 A(er) can be roughly estimated by setting 
the right hand side of Q2(J[1 equal to 1/2 of its maximal 
size, i.e 



Co<J~ \C 2 \(k(a) + A(a)) 2 



\C 3 \k d (a) 



hence 



A(a)«fc( ( x) 2 |C , 3|(2|C 2 |)- 1 =0( ( T) 



C. Estimating the decay time s,i 



(22) 



Let us first approximate U by U Sd , i.e. the time prop- 
agator for the time independent Dirac operator present 
at Sd- 

_3 

Using 1(131) 1 lUD an( i l|2*2*|l we have Sd oc es d 2 , hence 
Sd oc e=. The rigorous estimate, taking into account the 
time dependence of the external field, yields Sd oc £3 
([T^L Hence if the field stays overcritical for times S ^> 
es the probability of pair creation is one. Note that in the 
adiabatic case S — 0(1), and thus this is well satisfied. 

The distribution of the outgoing momenta will be dis- 
cussed below. 
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IV. e IN HEAVY ION COLLISIONS 

One way to create overcritical fields experimentally are 
heavy ion collisions. (There are other experiments, which 
might become more relevant for SPC E3)- 

bmce 

the fields of the nuclei are repulsive for positrons, it is 
the positron which scatters. For heavy ion collisions the 
adiabatic time scale on which the field increases is di- 
rectly determined by the relative speed with which the 
heavy ions approach each other and one computes that e 
is of order 10 _1 , however the time where the field remains 
overcritical is very small (see [12]). The time variation 
of the field (even for weak fields) produces also pairs (see 
e.g. [23] ) , which may become relevant when the field 
strength is close to criticality and where the time dura- 
tion of over-criticality is small (see below). In principle 
an estimate of the induced pairs is needed or an exper- 
imental measurement of this "background" : measuring 
once a system which is slightly under-critical and com- 
pare the rate of created pairs to a system which is slightly 
over-critical. 

Theoretical models indicate that the collision time can 
be enlarged if the nuclei form a composite nuclear system 
[2^|. This effect may be useful to increase the amount 
of spontaneously created pairs in U-U scattering experi- 
ments, though it seems not possible to enlarge the colli- 
sion time beyond the decay time of the previous chapter. 



We apply now P 1 - and replace in the integral U by Uq, 
which yields by iteration an error of smaller order. Ob- 
serving further that P ± $ = and <f> = D $ 

\\P ± U(S, 0)$|| 2 

= \\P ± - / U (S,s)(l - D s )U (s,0)$ds\\ 2 
6 Jo 

= \\ % -e=¥ / P^{l-D s )$dsf 
e Jo 

= \\-e^ [ P x (Ao-A s )<i>ds\\ 2 



This differs from estimates in the literature (see for ex- 
ample formula (20) in [24|), where the formula different 
from but corresponding to J3U|)(see e.g. formula 2.24 in 
|22|) is interpreted as a Breit Wigner form, leading to 
false estimates. We note that for the "short time anal- 
ysis" (which is S <C Sd) the use of Breit- Wigner form is 
not reasonable anyhow. 



A. SPC probability for Collision times <C Sd 

In the literature the SPC-probability has been com- 
puted for an overcritical field of very short life times like 
S <§ ••< -22]. 

We note that this probability can be easily estimated 
without any reference to the resonance, contrary to [4j. 
We also note, that for such short life times the adia- 
batically changing field can only vary by O(S). How- 
ever in the literature one considers also the case where 
the field changes — when it reaches over critical values — 
nonadiabatically to a very small overcritical stationary 
value of the size Aq + a with a«S. In this case only a 
very small part of the critical wave function $ will scatter 
in the positive continuous subspace. Let P^ be the pro- 
jection onto the subspace orthogonal to <E>. The probabil- 
ity for scattering may be estimated by \\P ± U(s, 0)$|| 2 . 
Now let Uo(s,0) — exp— then 

{U(S,0)-U o (S,0))$ 

= - / U {S,s)(l-D s )U{s,0)$ds . 
e Jo 



V. ENERGY SPECTRUM OF THE OUTGOING 
POSITRONS 



An interesting prediction is the shape of the momen- 
tum distribution of the created positron. It would be 
nice if the shape would be simply the resonance (|20|l as 
suggested e.g. in [22|. But that requires a somewhat dif- 
ferent situation than what has been discussed here. It 
would require an overcritical static field (adiabatic is not 
enough) of a life time S larger than the decay time s^. 
Then in fact the resonance H2U|) would stay more or less 
intact (see also 24] ). In an adiabatically changing field 
the resonance ()2(J|1 changes however with the field and it 
is highly unclear how. The resonance l|2U|l will be surely 
washed out and one can only expect a shape which is 
somewhat peaked around small momenta. 
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